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ABSTRACT 

A  power  transformation  regression  model  Is  considered  for  ex- 
ponentailly  distributed  time  to  failure  data  with  right  censoring, 
procedures  for  estimation  of  parameters  by  maximum  likelihood  and 
assissment  of  goodness  of  model  fit  are  described  and  illustrated. 
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1 .  INTRODUCTION 

Suppose  that  ve  have  n  individuals  or  components  which  are  subject 
to  'failure1  and  that  the  response  variable  of  interest  is  time  to 
failure.  Let  Tj*,  T2*,  ....  Tn*  be  independent  random  variables 
representing  the  times  to  failure  of  the  individuals.  We  suppose 
that  right  censoring  may  occur  because  of  the  need  for  early  termination 
of  the  investigation  and  let  Ti  ,  T2 , . . . ,  Tn  represent  the  recorded 
survival  times.  Defining  censoring  indicator  variables 


fl  if  T  *  i 
«|  *  1 


s  uncensored 


(1.0 


(0  if  T>  is  censored 


we  have 


We  let 


T.*  -  T.  if  u.  ■  1  ,  T.*  >  T.  if  co.  -  0 

1  1  1  1  1  i 


n 

n  -  /  u. 

U  i-1  1 


n  -  J  (J  -to  ) 
i-1  1 


(1.2) 


(1.3) 


denote  the  numbers  of  uncensored  and  censored  observations,  respectively. 


We  now  suppose  that  measurements  are  available  on  k  explanatory  variables 


xj  ,  x2,...,  xfc.  Setting  x'  -  (Xj  ,x2 , . . . txfc)  the  p.d.f.  and  survival 
function  of  T*  given  x  will  be  denoted  by  f(t;x)  and  S(t;x), 
respectively.  If  the  hazard  function  f (t;x)/S(t;x)  is  assumed  to  be 
a  constant  independent  of  t  for  any  given  x,  then  T*  has  the 
exponential  distribution  with  p.d.f. 


ri 


exp(-t/py),  t  >  0 


"  1  0  otherwise 

where  p^  denotes  the  expected  value  of  T*  given  x. 

Various  models  have  been  proposed  in  the  literature  to  represent  the 
dependence  of  p^  on  x.  Fiegl  and  Zelen  (1965)  consider  the  model  form 


(1.4) 


V  -  X(l*x’B) 


(1.5) 


I 


I*  • 

» 


r  ^ 


a 


vhile  Greenberg  et. al. (1974)  use  Che  form 


Ux  -  A/O  4  x'B) 


0.6) 


where  8*  -  (8.  ,8,, .  . .  ,8.  )  and  A  is  a  positive  constant.  Both  models 

*N.  12  K 

require  that  the  condition  x'B  <  1  must  be  imposed  to  ensure  that 
Wx  >  0.  An  alternative  model  which  does  not  require  a  constraint  to 
be  imposed  on  x*£  is 


Xexp(x'B) 


0.7) 


■  *'  i*' 

.  • .  .1 


This  model  arises  for  the  exponential  case  from  the  well-known  family 
of  proportional  hazard  regression  models  (see  e.g.Kay  (1977))  in  which 
an  assumed  underlying  hazard  function  is  adjusted  by  multiplicative 
exponential  factors  to  allow  for  the  effect  of  the  explanatory  variables. 
Prentice  0973)  also  discusses  the  use  of  censored  regression  models  for 
the  exponential  case. 

In  this  paper,  we  consider  the  power  transformation  model  given  by 


Px  -  X(1  4  Bx'g)1 


♦  6x'B)l/6 


0.8) 


We  shall  refer  to  $  as  the  power  parameter.  It  is  seen  that  when 
6*1,  the  model  corresponds  to  the  usual  additive  regression  model 
for  the  mean  given  by  (1.5),  while  if  6  *  -1  the  reciprocal  model 
(1.6)  is  obtained  after  appropriate  reparameterisation.  When  6 -4  0, 
the  model  given  by  (1.7)  is  obtained.  For  later  work,  it  is  convenient 
to  write  the  power  transformation  model  in  the  form 


(1  ♦  6x’B  )*/6 


0.9) 


where  *;  -  0,x’),  8;  -  (B^.B^,. 


(A6-l)/6,  B..  -  X6B.,  j-1 . k 

J  J 


0.10) 


In  general,  the  power  parameter  6  as  well  as  the  coefficient  vector  £ 
will  have  to  be  estimated  from  the  data.  In  section  2,  results  are  given 
for  maximum  likelihood  estimation  and  it  is  shown  how  the  estimates  can 
be  obtained  numerically  using  the  statistical  package  GLIM.  The  problem 
of  assessment  of  the  goodness  of  fit  of  the  model  is  discussed  briefly  in 
section  3  and  in  section  4  a  numerical  example  is  used  to  illustrate  the 
procedure. 


*  -*  -  -  ■  ■  -  •  -  »  j  - 


,V 


2. 


MAXIMUM  LIKELIHOOD  ESTIMATION 


We  shall  assume  that  the  stochastic  model  underlying  the  censoring 
mechanism  is  unknown.  The  likelihood  given  the  data 

where  x?  -  (xii»xi2»* ,,,xik>*  is 


nr  1  1-w. 

L  -  n<f(t.  ;  x.  )}  1  {S(t.  ;  x.)f 

i.jl  1  1  x  1J 

”  n  {uf1  exp(-t./y  )} 
i-ll  ii  1  *i  J 


where  U.  *  (1  ♦  ) 


»/6 


and  x*.  -  <l,x\).  If  we  put 

1  —  1 


-i 

x. 

~i 


1*1 | 2| • • • |D 


the  log  likelihood  is 


log  L  -  J  (u.logK  -  D.)  -  j  “.logt. 
i-1  1  1  1  i-1  1  1 


The  likelihood  equations  are  obtained  by  setting 

n  <Vy.)  3Vi  _ 

- w. 


n  (0..-V.)  3Vi  .  o  " 

l  ^  .1 


i-1  Vi 


j  »  0, 1 , . .  .k 


i-1  yi 


*3 


Since 


the  maximum  likelihood  estimates  are  found  as  the  solution  of  the 
k  ♦  2  equations 


«  *  /Vi\« 

j  x  ^ui”^iHt"")  . ^ 

i-1  ij  1  1  VV 

5  .  *  .  it  yi\*  .  /Pi  M  . 


I, 


where 


Ui  -  t.d  + 


Routine  calculation  gives  the  second  order  derivatives  of  the  log- 
likelihood  as 


SfogL  -  -  l  Jt.  x...  {  (6+  1)  u.  -  6ux  }  -  -  a.., 

31.  SSL,  i-j  1  >•  1  XJ  JJ 


32lo: 

"5K 


"  “  ij,  xij®i{u>ryi)cV,) "  Vi}  *  _bj 


^  "  ‘  i2.^[wizi+  (vyi){*i(3_®i)  •  ~V“  lo»0i}] 


“  -c 


for  j ,3 '  ■  0, 1 , . . . ,k,  where 


•i"(r)  •  ‘tm  i 


Setting 


t-  «ajj.»  . 


B  -  ((6.))  ,  C  -  ((c)) 

J 


A  A  A 

where  a..,  denotes  a..,  evaluated  at  8  «  8  and  6  =  c»  etc,  the 

jj  jj  . 

estimated  asymptotic  covariance  matrix  of  the  estimators  is 


V 

V 

A 

B  ‘ 

V  - 

~n 

~li 

. 

V1 

V 

B’ 

C 

.  "”12 

~22  J 

*** 

l-l 


The  matrices  in  the  partitioned  form  are  given  by 


\  -  bl*  h'm  -  r*  V 1  s' a-1 


(2.11) 


l„-  -Vuc-i'h'v1  >  -  (z-z'b'v 


W  ■>  V  '.J.  '.J 


The  calculation  of  the  estimates  and  statistics  used  for  making 
inferences  about  6  and  £„  is  performed  straightforwardly  using  the 
statistical  package  GLIM  which  has  been  developed  for  fitting  generalised 
linear  models  by  maximum  likelihood  for  error  distributions  in  the 
exponential  family.  The  method  is  similiar  to  that  described  by  Aitken 
and  Clayton  (1980)  for  fitting  proportional  hazard  regression  models  for 
censored  survival  data. 

In  (2.3),  the  first  term  represents  the  kernel  of  a  log- likelihood 
treating  Wj ,  W2,...,Wn  as  independent  Poisson  variables  with  E(VL)  *y^, 
and  the  second  term  is  independent  of  B  and  6.  From  (2.2)  we  have 


{(t./y.)6  -  0/6  -  x*.  6  (2.12) 

X  1  * 

A 

If  6  is  known,  the  value  6  (6)  which  maximises  the  log-likelihood, 

now  denoted  by  logKg  (6),  6)  ,  can  be  found  fitting  a  Poisson  regression 

*  6 
model  with  the  user  defined  link  function  g. (V.)  » {(t./y.)  -  l}/6. 

’J  X  IX 

In  the  usual  case  when  6  is  unknown,  the  maximum  likelihood  estimates 

A  A 

£  and  £  can  be  found  using  the  following  iterative  procedure  which 
is  based  on  a  general  method  suggested  by  Pregibon  (1980)  for  selecting 
the  link  function  in  generalised  linear  models. 


Suppose  that 


5t  denotes  an  initial  approximation  to  5.  Using  the 


first  order  Taylor  series  expansion  about  the  value  6  ■  6j  gives  the 
approximation  a  a  a 

I,)®1-  1  a  r  K  <ti/u,)6‘  log(t./y.)  +  i  -  (t./y.)6> 

-  ♦  (6-6^1  — — i  - ^  } 


*6(V  * 


(ti/y.3 


& 


(2.13) 


An  initial  fit  is  made  using  the  model  g*  (U^)  ■  to  yield  the 

estimate  |*l)  and  fitted  values  fif*5  from  which  we  form 


r-  (2-,4) 


A  fit  of  the  model  g£  (M*  )  *  x'.f*  ♦  *•  (6  -  6)  then  gives  the  second 

a  a(i)  *x  1 

approximation  62,  where  ^  ,  is  treated  as  an  additional  explanatory 
variable.  The  procedure  is  then  repeated  until  satisfactory  convergence 
is  obtained. 


.•  V  VVV  '/'/'j*** 


If  a  confidence  interval  for  5  is  required,  we  use  the  approximate 
distribution  result 

2  log  L  {|M6) ,  6}  -  2  log  U  |^(6) ,  8  X? 
giving 

P  [logL{|^(6).6}>logL(|^(S),S)-ixl2<,“a)]  *  (2.15) 

where  X^O  ~ci)  denotes  the  upper  lOOoft  point  of  the  x^  distribution. 

Thus  if  6^  and  6^  are  the  two  values  of  6  for  which 

logL{B#(6),6}  -  logL{^($),$)  -  ixf C 1  “ a) 

where  6. <  C  ,  then  the  interval  (6t,6  )  provides  an  approximate 
L  U  L  U 

100(1-  a) Z  confidence  interval  for  6.  The  calculation  of  this  interval 

A 

will  require  values  of  logL{8*(5),6}  to  be  determined  for  a  grid  of 
values  of  6. 

A  A 

The  covariance  matrices  for  £  and  6  can  also  be  found  straightforwardly 

using  (2.11)  since  the  GLIM  output  will  provide  the  value  of  A-1, 

a  a  ~ 

the  covariance  matrix  for  when  6  is  treated  as  the  true  fixed  value 
of  6.  Finally,  we  may  wish  to  test  the  hypothesis  that  a  specified 
subset  of  the  explanatory  variables  have  no  effect  on  the  expected  time 
to  failure.  This  can  be  done  using  a  standard  likelihood  ratio  test  based 
on  the  difference  between  the  maximised  log-likelihoods  under  the  full  and 
reduced  models,  this  being  the  difference  between  the  corresponding 
deviances  in  the  GLIM  output. 

3.  RESIDUAL  PLOTS 

A  graphical  procedure  based  on  suitably  defined  observed  residuals  provides 
a  useful  preliminary  assessment  of  adequacy  of  fit  of  the  power  transformation 

.j/1 

model.  If  the  model  is  correct,  the  random  variable  R  ■  T*(l+6x’B  ) 
has  the  standard  exponential  distribution.  Thus  if  there  is  no  censoring, 
we  define  observed  residuals  by 


If  the  model  is  correct,  the  {r^}  will  approximately  have  the 

properties  expected  of  a  random  sample  of  n  observations  from  the 

standard  exponential  distribution.  A  plot  of  the  ith  ordered  residual 

n 

r,.N  against  the  ith  smallest  exponential  score  e.  *  ]  j  should 
M  *•"  n-i.i 

then  give  an  approximate  straight  line  relation  with  unit  slope  and  zero 
intercept. 


To  deal  with  censoring,  we  note  that 


E{R|  R>  t  (1  +  fix’g  )"  1/6}  -  1  +  t(1  +  fix '8  )' 


(3.2) 


Replacing  the  censored  value  of  R  by  the  maximum  likelihood  estimate 
of  it6  expected  value,  we  define  modified  residuals  by 


{a  A  _i/(5 

“  V' 

I  *t.(l  ♦  fix'  .6  )-1^  “  v° 

1  ~*1~* 


(3.3) 


A  plot  of  r. 


against  e.  should  again  give  an  approproximate 
i,n 


straight  line  relation  if  the  assumed  power  transformation  model  is 


correct. 


4.  AN  ILLUSTRATIVE  EXAMPLE 

To  illustrate  the  use  of  the  power  transformation  model  we  use  the 
data  in  table  1  which  were  taken  from  Prentice  (1973).  The  data 
are  for  groups  of  advanced  lung  cancer  patients,  the  groups  being 

defined  by  two  factors,  type  of  chemotherapeutic  agent  at  two  levels 
standard  and  test,  and  tumour  cell  type  at  four  levels,  squamous,  small, 
adeno  and  large.  In  the  original  study,  four  explanatory  variables  were 
measured  and  the  analysis  based  on  model  (1.7)  indicated  that  only  one 
variable,  general  medical  status  on  a  scale  10, 20,..., 90  had  a  real 
effect.  Ve  have  only  included  this  variable  in  our  analysis.  The 
survival  times  are  in  days  and  right  censored  observations  are  shown 
with  an  asterisk. 

The  GLIM  fitting  procedure  described  in  section  2  was  used  to  fit  the 

model.  .» 

M,(j)  -  [1  +  6(8  (j)  +x8  (j)}]  /0,  j-1,2 . 8  (4. 


If  the  model  is  correct,  the  (r-J  will  approximately  have  the 

properties  expected  of  a  random  sample  of  n  observations  from  the 

standard  exponential  distribution.  A  plot  of  the  ith  ordered  residual 

n 

r,^  against  the  ith  smallest  exponential  score  n  “  I  J  should 

*  n-i+I 

then  give  an  approximate  straight  line  relation  with  unit  slope  and  zero 
intercept. 


To  deal  with  censoring,  we  note  that 


E{R|R>to  +  6x;^)‘l/6}  -  i4t(i*g»;ex)~vg 


(3.2) 


Replacing  the  censored  value  of  R  by  the  maximum  likelihood  estimate 
of  it6  expected  value,  we  define  modified  residuals  by 
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r? 

1 


> 


A  A  -i/fi 

’  A  if  “i-i 

i  >'1/S  if  V° 

1  ~*1~* 


(3.3) 


A  plot  of  r,..  against  e.  should  again  give  an  appropr^xi  te 
UJ  i«o 

straight  line  relation  if  the  assumed  power  transformation  model  X6 
correct . 


4.  AN  ILLUSTRATIVE  EXAMPLE 


To  illustrate  the  use  of  the  power  transformation  model  we  use  the 
data  in  table  1  which  were  taken  from  Prentice  (1973) .  The  data 
are  for  groups  of  advanced  lung  cancer  patients,  the  groups  being 

defined  by  two  factors,  type  of  chemotherapeutic  agent  at  two  levels 
standard  and  test,  and  tumour  cell  type  at  four  levels,  squamous,  small, 
adeno  and  large.  In  the  original  study,  four  explanatory  variables  were 
measured  and  the  analysis  based  on  model  (1.7)  indicated  that  only  one 
variable,  general  medical  status  on  a  scale  JO, 20,..., 90  had  a  real 
effect.  We  have  only  included  this  variable  in  our  analysis.  The 
survival  times  are  in  days  and  right  censored  observations  are  shown 
with  an  asterisk. 

The  GLIM  fitting  procedure  described  in  section  2  was  used  to  fit  the 

modei.  ... 

lix(j)  -  tl  *6{B^(j)  ♦*^1(j))J  j*l  ,2, . . .,  8  (4.1) 


TABLE  1 

_ i 

• 

j 

Data 

for 

lung  cancer 

patients: 

t  - 

days 

of  survival. 

X  * 

general 

Y.  1 
..  *•] 

medical 

status. 

G, 

*  standard 

squamous. 

G2 

■  standard 

small , 

G3  - 

standard  adeno. 

Gw 

■  standard 

large,  G 

s  * 

test. 

squamous. 

Gt  * 

test,  small,  G 

7  * 

test. 

adeno 

»  gb 

“  test, 

large. 

• 

Gi 

g2 

G 

3 

Gw 

c5 

g£ 

g7 

Ge 

t 

X 

t 

X 

t 

X 

t 

X 

t 

X 

t 

X 

t 

X 

t 

X 

72 

60 

30 

60 

8 

20 

177 

50 

999 

90 

25 

30 

24 

40 

52 

60 

-- 

4U 

70 

384 

60 

92 

70 

162 

80 

112 

80 

103* 

70 

18 

40 

164 

70 

• 

228 

60 

4 

40 

35 

40 

216 

50 

87* 

80 

21 

20 

83* 

90 

19 

30 

126 

60 

54 

80 

117 

80 

553 

70 

231* 

50 

13 

30 

31 

80 

53 

60 

1  18 

70 

13 

60 

132 

80 

278 

60 

242 

50 

87 

60 

51 

60 

15 

30 

10 

20 

123* 

40 

12 

50 

12 

40 

991 

70 

2 

40 

90 

60 

43 

60 

82 

40 

97* 

60 

162 

80 

260 

80 

111 

70 

20 

30 

52 

60 

340 

80 

- 

- 

1 10 

80 

153 

60 

3 

30 

200 

80 

1 

20 

7 

20 

73 

60 

133 

70 

• 

j 

314 

50 

59 

30 

95 

80 

156 

70 

587 

60 

24 

60 

8 

50 

1 1 1 

60 

100* 

70 

117 

80 

182* 

90 

389 

90 

99 

70 

36 

70 

231 

70 

42 

60 

16 

30 

143 

90 

33 

30 

8 

80 

48 

10 

378 

80 

8 

40 

151 

50 

105 

80 

25 

20 

99 

80 

7 

40 

49 

30 

. '  :i 

144 

30 

22 

60 

103 

80 

357 

70 

61 

70 

140 

70 

* ' 

-  *'-H 

25* 

80 

56 

80 

250 

70 

467 

90 

25 

70 

186 

90 

1. 

1 1 

70 

21 

40 

100 

60 

201 

80 

95 

70 

84 

80 

18 

20 

1 

50 

80 

50 

19 

50 

•  i 

139 

80 

30 

70 

51 

30 

45 

40 

- 

20 

30 

44 

60 

29 

40 

80 

40 

-  ■*-  ■, 

31 

70 

283 

90 

— 

52 

70 

15 

50 

i 

'  .-r* 

287 

60 

/■ 

18 

30 

51 

60 

Yv*- 

122 

80 

,\ 

27 

60 

m 

54 

70 

■Y 

—  — ! 

7 

50 

63 

50 

392 

40 

10 

40 

I 

where  Wx(j)  denotes  the  expected  survival  time  for  a  patient  in  group  j 

having  medical  status  value  x.  The  use  of  the  initialisation  value  -\_ 

a  aaa 

*  0  led  to  the  sequence  62  ■  0.359,  63  *  0.432,  64  ■  0.429  of  ‘ 

A  A  - 

approximations  to  6  and  6  was  therefore  taken  as  0.43  correct  to  V 

.  a 

two  decimal  places.  In  table  2,  values  of  V>x(j)  are  given  for  the 

eight  groups  and  the  nine  medical  status  levels  used  for  classifying 

the  patients.  The  maximum  likelihood  estimates  of  B#0(j)  and  0^(j) 

are  also  shown  together  with  their  estimated  standard  errors.  The  results  l~_ 

show  that  at  the  higher  levels  of  x,  the  expected  survival  times  are 


I 


3 


highest  for  group  G5,  and  that  the  differences  among  the  groups 
G3,  G6  and  G7  are  relatively  small.  These  results  support  the  general 
finding  made  by  Prentice  and  it  should  be  noted  that  for  these  data, 
the  use  of  the  power  transformation  model  leads  to  only  a  moderate 
improvement  in  fit  compared  with  the  model  given  by  (1.7).  The  value 
of  2  log  L{  (j5^ (0.43) ,  0.43)  -  2  log  (0)  ,0)  is  1.41  which  when 
referred  to  the  Xj  distribution  is  just  significant  at  the  252  level. 

It  is  of  course  instructive  to  use  the  general  power  transformation  model 
to  assess  the  adequacy  of  fit  of  the  commonly  used  model  given  by  (1.7). 

TABLE  2 

Maximum  likelihood  estimates  of  regression  coefficients  and  of  expected 
survival  times . 


X 

g2 

g3 

g5 

g6 

g7 

G6 

10 

27.4 

43.0 

1.4 

49.2 

1.9 

9.3 

23.8 

2.5 

20 

43.6 

52.1 

5.6 

67.4 

13.8 

15.5 

30.1 

10.3 

30 

64.1 

62.2 

13.0 

88.8 

39.4 

23.6 

37.1 

24.8 

40 

89.2 

73.4 

24.2 

113.7 

81.2 

33.7 

45.1 

46.8 

50 

119.2 

85.6 

39.5 

142.1 

140.7 

46.0 

53.9 

77. 1 

60 

154.2 

98.9 

59.2 

174.3 

219.4 

60.4 

63.6 

116.2 

70 

194.5 

113.3 

83.5 

210.2 

318.6 

77.1 

74.3 

164.7 

80 

240.1 

128.8 

112.7 

250.0 

439.4 

96.1 

85.9 

223.0 

90 

291.4 

145.5 

146.9 

293.7 

582.7 

117.6 

98.5 

291.7 

k0<v 

5.21 

8.38 

-1.74 

8.30 

-3.36 

2.23 

5.81 

-1 .80 

SE{B#0(j)}  6.76 

4.83 

2.61 

13.75 

3.23 

2.68 

2.86 

4.48 

0.213 

0.101 

0.214 

0.179 

0.411 

0.150 

0.096 

0.291 

SE{6^1(j)}  0.140 

0.090 

0.072 

0.204 

0.999 

0.067 

0.055 

0.110 

5.  EXTENSIONS 

The  regression  model  given  by  (1.8)  for  the  mean  or  equivalently  the 
reciprocal  of  the  hazard  function  of  the  exponential  distribution  is 
more  flexible  than  other  models  that  have  been  previously  used.  It 
is  therefore  likely  that  except  for  very  large  data  sets,  the  model  will 
provide  a  satisfactory  fit  if  the  underlying  distribution  is  exponential, 
even  if  the  model's  link  function  is  not  correct. 

For  distributions  other  than  the  exponential  the  power  transformation 
link  function  may  be  used  to  model  the  dependence  of  the  hazard  function 


-•  .>  .•* , 


on  explanatory  variables.  Thus  if  h(t;  x)  denotes  the  hazard 
function  for  tine  to  failure  given  :x,  we  may  use  the  proportional 
hazards  regression  model 


h(t;x)  «  hc(t;0)(l  +  6x’B)1'6  (5.1) 

where  h  (t;8)  represents  che  hazard  function  for  an  individual  with 
x  *  0  which  is  assumed  to  depend  on  a  vector  6  of  unknown  parameters. 


A  straightforward  argument  shows  that  the  log-likelihood  under  the  model 
given  by  (5. 1)  is 


log  L  (8,6,6) 


n  n 

I  («iiogul-yi)+I  wilog{h6(t.;e)/Ho(ti;e)} 
i»l  i-1  A  '  ‘  ' 


(5.2) 


where 


H  (t;6)  ■  f  h  (u;8)  du 

®  I  O  M 


(5.3) 


is  the  cumulative  hazard  function  and 

p.  -  (I  +6x!6)1/5H  (t;6) , 

1  fwl  Q  rw  9 


i»I, . . .  ,n 


(5.  A) 


Suppose  now  that  the  underlying  distribution  depends  on  two  parameters, 
X  and  6  say,  and  that 


Ho(t;X,6)  -  Xn(t;6) 


(5.5) 


where  X  >  0  and  n(t;0)  2  0Vt>0.  Examples  of  distributions  having 
this  property  are  the  Weibull  with  Ho(t;X,8)  «  Xt  and  the  Gompertz 
with  Ho(t;X,6)  ■  (X/8)exp(0t) ,  t>0.  If  (5.5)  holds,  we  may  write 


P.  -  n(t.;8)(l  +  6x\B  )J/5 

i  i 


where 


j 


1 , . . . ,  k 


(5.6) 


(5.7) 


For  8  known,  the  ML  estimate  8  (8)  can  be  found  treating  the  {u>.} 

1 

as  observations  on  independent  Poisson  regression  variables,  the  link 
function  being  . 

.<V  (5.8) 

If  8  is  unknown,  a  search  over  a  grid  of  values  of  6  would  be  needed 
to  find  the  global  maximum  of  the  log-likelihood. 
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